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ABSTRACT 

m 
> 

^ I Emission line galaxies are an invaluable tool for our understanding of the evo- 

lution of galaxies in the Universe. Imaging of deep fields with narrow-band filters 
allows not only the selection of these objects, but also to infer the line flux and 
the equivalent width of the emission line with some assumptions. The narrow- 
band filter technique provides homogeneous samples of galaxies in small comoving 



O 



Qh' volumes in the sky. We present an analysis of the selection of emission-line galax- 

O ' ies using narrow-band filters. Different methods of observation are considered: 

^ • broad-band - narrow-band filters and two broad-band and one narrow-band fil- 

■ ters. We study also the effect of several lines entering simultaneously inside the 
filters (this is the case of Ha) . In each case the equations to obtain the equivalent 

^ ' width and line flux from the photometry are obtained. Candidates to emission- 

■ line objects are selected by their color excess in a magnitude-color diagram. For 
different narrow-band filters, we compute the mean colors of stars and galaxies, 
showing that, apart from galaxies, some types of stars could be selected with 
certain filter sets. We show how to compute the standard deviation of the colors 
of the objects even in the usual case when there are not enough objects to deter- 
mine the standard deviation from the data. We present also helpful equations to 
compute the narrow-band and the broad-band exposure times in order to obtain 
minimum dispersion in the ratio of fluxes of both bands with minimum total 
exposure time. 



Subject headings: techniques: photometric — stars: emission-line. Be — galaxies: 
fundamental parameters 
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Introduction 



Emission line galaxies (hereafter ELGs) are an invaluable tool for our understanding of 
the evolution of galaxies in the Universe. Faint galaxies can be hard to confirm spectroscop- 
ically, while the ELGs are generally easy to identify. Furthermore, the emission lines are 
produced within regions related either with star formation or with the active galactic nuclei 
(AGN) phenomena. 

Before the advent of CCD detectors, surveys for emission- line_galaxies used photo 



graphic plates with obiective prisms in Schmidt telescopes (e.g. 


Smith 


1975; 


Smith et al. 


1976 


MacAlDine et al 


Il977: Wasilewski 1983: Pesch & Sanduleak 


1983; 


Salzer et al. 


1989; 


Zamorano et al. 


1994. 


19961). 



Several low redshift surveys have made use of the objective-prism technique with CCD 
detectors to circumven t the limit on wave l ength and quantum efficiency imposed by pho - 



tographic plates (e.g., ISabbey et al.l l200ll ; iGronwall et al.l l2004l : iBongiovanni et al.l l2005l ) . 



However the objective prism, as any system of slitless spectroscopy, blends the flux of the 
object with that of the dispersed sky, reducing the sensitivity of the survey for the objects 
with less contrast between the continuum and the line (that is, with small equivalent width) . 
This effect, together with relatively small aperture of the telescopes used, limits these surveys 
to detect bright sources in a near volume of the space. 



A more efficient approach to detect ELGs is to use narrow-band iniaging (IDjorgovski et al 



19851 ; iLowenthal et al.l Il99ll ; iMacchetto et al.l Il993l ; iThompson et al.l Il995l ) . This technique 
reduces significantly the contribution of the sky brightness, as it is admitted in a small range 
of wavelengths. The sky background, that is the most significant limitation in the detection 
of objects in deep broad band images, is greatly reduced. A small wavelength range increases 
also the contrast between the emission line and th e continuum. Narrow-ba nd imaging has 
been used also to detect both galac t ic (lRauchlll999l) as we ll as extragalactic (lOkamura et al. 



2OO2I ; ICastro- Rodriguez et al.ll2003l ; lArnaboldi et al.ll2003l ) planetary nebulae. 



In order to use the narrow pass-band filters with maximum efficiency, regions of the 
night sky spectrum with a minimum background are selected. In the optical wavelengths, 
the windows in the Meinel OH bands around 8200A and 9200A have been used to center 
narrow-band filters. Most of the interest has gone to the region around 8200A, where the 
CCD detectors still have a good quantum efficiency. 



Similar to the narrow-band imaging is the use of a tunable ffiter (TF, iBland-Hawthorn fc Jones 



19981 ). A tunable ffiter works like a Fabry-Perot interferometer modified to produce a wide 
region of interference and to cover a range of wavelengths. 
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Compared with narrow-band imaging, a TF provides a pseudo-spectrum, when the 
spectral range is explored with small wavelength steps. As a disadvantage, the TF requires 
a dedicated instrument, whereas narrow-band imaging demands only the narrow-band filter 
and can be carried out in any imaging facility. Both narrow-band imaging and TFs need 
very long exposure times to achieve the necessary signal-noise ratio. 

The search for emission-line galaxies using a given narrow-band filter is open to different 
lines at different redshifts. Lines whose redshifted wavelength fall inside the narrow-band 
filter range could be detected. When the emission-line is bright enough the galaxy will be 
selected. The more prominent lines that can be detected in the 8200A and 9200A atmospheric 
windows are HaA6563 (at redshift z ~ 0.24 and z ~ 0.4), [Olll] AA4959,5007 (z ^ 0.6 and 
z ~ 0.8), [Oll]A3727 {z ~ 1.2 and z ^ 1.5) and LyaA1216 (z ^ 5.7 and z ~ 6.5). 



The Ha line is an exc ellent tracer of the Star Formation Rate (SFR, lKennicuttlll998 



Chariot fc Longhettill200ll ). provided that the ionizing flux comes from young stars and not 
from non-thermal activity. Other commonly used SFR tracers are the far infrared (FIR) and 
the UV. The three tracers share sensitivity to the parameters of the Initial Mass Function. 
Ha SFR is affected by obscuration but is not very sensitive to metallicity. The FIR is not 
affected by obscuration, but there is uncertainty in how the total (8-1000/im) luminosity is 
computed from monochromatic measurements (e.g. Spitzer's 24/im) and the contribution of 
different non-SFR related components (cirrus, ho t dust) to the FIR. Finally the UV is heavily 
affected by obscuration. As shown by iBelll (120031 ) . obscuration corrected Ha is consistent, to 
a factor of 2, with the summed SFRs estimated using UV and FIR(8-1000/im). Consequently, 
Ha observations of galaxy samples with UV and FIR data provides an invaluable tool to 
understand the evolution of the SFR and the role of obscuration in the determination of 
global SFR for galaxies. Some works have us e d narrow-bari d imag i ng to select Ha candi dates 
in the 8200A wind ow: IPascual et all J200lh : iFuiita et all J2OO3I ): lUmeda et aP tooi ) and 
Aiiki et al.l J2006bl l 



Given that small amounts of dust reddening can completely extinguish Lya-emission, 
Ha is the strongest emission line in many starburst galaxies. The [Olll] AA4959,5007 doublet 
is in many cases the other strongest emission line. In fact, the [Olll] AA4959,5007/Ha ratio 
is 0.6 in local star-forming galaxies (IMouhcine et al.ll2005l ). This line is not used to trac e 
the star formation but to trace the power of AGN (see, for example iKauffmann et al.ll2003l ). 
Thus, [Olll] AA4959,5007 emitters can be a tool to unveil AGN at different redshifts. 

[On] A3727 i s also intense when compared with Ha, with a ratio [Oil] A372 7/Ha = 0.45 



(IKeiinicuttl Il992[ ) . However, the ratio depends on luminosity and metallicity (jJansen et al. 
2OOII ) and varies from sample to sample. This line is used as a star-formation indicator at 
redshifts where Ha is outside the visible range (which happens at z > 0.4). SFR based 
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on [Oll]A3727 comes from the fact that there is a good correlation between [Oll]A3727 



cahbrations of SFR based on Ha are used (see , e.g. 



Gallagher et al. 


1989; 


Kennicutt 


1998 



Rosa-Gonzalez et al.ll2002l : iKewley et al.ll2004j ) but, in general, the correlation between S FR 
and [Oll]A3727 is affected intensely by extinction and metallicity (IMouhcine et al.l 120061 ). 



In the assumption of case B recombination, two-thirds of the Lyman continuum pho- 
tons are reprocessed as Lya. Therefor, young starbursts should be easily detected by 



their Lya emissiori. How ever, observations of nearby starbursts (e.g. iHartmann et al.lll988 



Calzetti fc Kinneylll992l ) revealed in most starburst galaxies a much weaker Lya than pre- 
dicted by simple models of galaxy formation. To explain this low flux, I Chariot fc Fall! (119931 ) 
sugg ested that Lya strong em ission requires either a low abundance of dust or an AGN. How- 



ever, 



Giavalisco et al.l (Il996l ) demonstrated that pure extinction could not explain the weak 
Lya fluxes. Results suggest that Lya pho tons are decoupled from the continuum radiation 
resulting in a increased sensitivity to dust iNeufeldl (Il990l ). 



Lya systems a.t z > 5 are a very useful probe of the properties of very young star-forming 
systems at the epoch after reionization. Unfortunately , the data obtained to the date (e.g. 



Kodaira et al 



2005 



2003 



Ajiki et al 



IHu et al 



2006a 



2004 



Rhoads et al.ll2004l : IWestra et al.ll2005l : iTaniguchi et al. 



Shimasaku et al.ll2006l ) are not sufficient to place strong constraints 



on the fundamental properties. The Lya line is not used commonly as a SFR tracer. Res- 
onant scattering of Lya photons by ne utral atomic hydrogen affects their relation with the 
SFR in a galaxy Jcharlot fc Falllll993h . 



Narrow-band imaging produces approximately volume- limited samples, since the narrow- 
observed bands correspond to small windows in redshift space. The objects are selected with 
a well defined limit in equivalent width and the line fiux can be transformed into luminosity 
with some simple assumptions. Narrow-band imaging thus provide line luminosities for a 
volume-limited sample of emission line galaxies. In the case of detecting lines used as star 
formation tracers, the sample would be directly SFR-selected, except for the AGN contribu- 
tion. The problem with this approach is that stars contaminating the sample or contributions 
from different emission lines cannot be separated only wit h narrow-band imaging. Additio nal 
assumptions (on the luminosity functions of ELGs, e. g. iJones &: Bland-Hawthorrul200ll ) or 



additional data (see, e.g., the color-color diagrams of iFuiita et aDl2003h are needed to com- 
plete the scenario. 

The aim of this work is to contribute to several aspects of the narrow-band imaging in 
order to optimize the results of the observations. Section [2] describes the characterization of 
the filters used through this work. 
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In Section [3] we deal with the different sources of contamination, either galactic or ex- 
tragalactic that can be found in an narrow-band survey and different methods used to detect 
them. Section H] describes how to infer the line and continuum fluxes (and the correspond- 
ing equivalent width) from the measured magnitudes in a variety of scenarios, including 
more than one broad-band filter and more than one line inside the narrow-band filter. We 
study the method to select the candidates, based solely on broad and narrow-band colors 
in Sect. El We study in detail also the comoving volume covered by a narrow-band selected 
sample. Numeric examples with three filter set, including broad and narrow-band filters, are 
shown in Section [6l Although we have based our examples in two surveys carried out using 
filters at the 8200A and 9200A windows, the analysis and methods described in this work 
are valid for a wide variety of narrow-band imaging surveys, including those extended into 
the JR. Finally, in Sect. [3 we present helpful equations to compute the narrow-band and 
the broad-band exposure times in order to obtain minimum dispersion in the ratio of fiuxes 
of both bands with minimum total exposure time. Conclusions are exposed in Sect. [HI 



Characterization of the filters 



In general, the filters used in CCD photometry are characterized by their transmittance 
curve T(A). The fiducial transmittances are usually measured by the observatories with col- 
limated beams and at room temperature. The total transmittance includes the contribution 
of the quantum efficiency (QE) of the detector and the transmittances of the optical systems 
in the light path. The working temperature, usually lower than the room temperature, pro- 
duces a blue shift of the transmittance of the narrow band filters. A converging light beam 
passing through the filter produces a second blue shift, as well as a decrease of the overall 
transmittance and a broadening of the transmittance profile. These effects depend on the 
focal ratio of the beam passing through the filter. 

Except in Section [5], where absolute source fiux are needed for statistics, we use always 
the transmittance normalized to peak unity (T) instead of the absolute transmittance (T). 
Both are related by T = TpeakT, where Tpeak is the maximum value of the transmittance. 

The mean wavelength of the band, An, is computed in agreement wit h the conventional 
definition (see for example iFukugita et al.lll995l : iFiorucci fc Munarill2003r ). 



Ar 



/Ar(A) dA 



A is the width of the filter computed from T. It agrees with the definition of the equivalent 
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width of the band transmittance profile Wq of iFiorucci fc Munaril (120031 ): 



Wo = A = j r(A) dA 
We use A instead of Wq and reserve W for the equivalent width of the emission lines. 



(2) 



/i^ (Eq. 7 of iFiorucci fc Munaril |2003| ) is the second-order central momentum of the 



transmittance profile: 



fx 



/(A-Ao)^T(A)dA 
/r(A)dA 



(3) 



The fiux calibration of the narrow-band images can be achieved using spectrophoto- 
metric standard stars, by obtaining synthetic magnitudes using the known SED of the spec- 
trophotometric stars. 

There are different choices for the election of the zero point of the narrow-band filter 
magnitude. The particular election does not affect the selection of objects or the measured 
fluxes. Different zero points in the bands only introduce an offset in the broad-narrow color. 

One approach is to use AB magnitudes (IOkdll974l ) or ST mag nitudes (J. Walsh 1995)0 
that provide a flxed physical magnitude scale. In other approximations, a broad-band fllter 
is chosen as reference and the zero point of the narrow band is obtained by flxing a zero color 
for a particular SED. We adopt, unless noted otherwise, the normalization of color zero for 
the flat fx spectrum. This is equivalent to make equal the zero points of both bands. 



3. Contaminating sources 

The redshift of each source entering the narrow-band fllter is unknown a priori. Sources 
with strong emission lines (particularly, but not only. Ha, [Olll] AA4959,5007, [Oll]A3727 
and Lya) produce a narrow-band flux excess. Different emission lines at different redshifts 
cannot be distinguished using only the narrow-band excess emission. Additionally, other 
sources (either galactic or extragalactic) without emission lines can exhibit a narrow-band 
flux excess when their SEDs passes through the fllter set used for ELG selection. Galaxies and 
stars with abundant absorption features in the atmospheric windows studied can potentially 
appear in narrow-band surveys as emitting objects. This aspect has to be studied for each 
fllter set particularly (in Section [6], examples with different fllters are shown). 



Optical and UV spectrophotometric standard stars ^http : //www . eso . org/ observing/ stemdards/ spectra 
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In general, stars show colors similar to those of the black body of their effective tempera- 
ture. For the late spectral types, the presence of molecular bands change the trend and large 
colors can be observed. In the specific wavelength range of the two atmospheric windows, 
the absorption produced by some molecular species present in the atmosphere of cold stars 
(TiO, VO) dominates. For late K stars, broad-narrow colors diverge from the black body. 
White dwarfs with broad Hydrogen absorption lines would produce also a narrow-band flux 
depression. 

It is important to note that depending on the filters used, particularly on the relative 
position of their Aq, molecular bands of cold stars could produce an increase of the color. 
Those stars with a large color excess produced by absorption lines can be misclassified 
as ELGs. As the ratio of sta rs to galaxies decreases with increasing magnitude (see, for 



example, iKaraali et al.l 120031 ) . the number of contaminating stars do not dominate the 



candidate selection at faint magnitudes. In any case, particular care has to be taken in the 
design or election of the filter set in order to avoid this source of contamination. 

Some types of stellar objects show emission lines in the wavelength interval where the 
filters are defined. Novae, symbiotic stars, cataclysmic variables (in general, low mass stars 
with either coronal activity or binaries) can show emission in the line Pa9 (9229 A). After 
a nova explosion there is emission in the multiplet Ol near 9261-9266 A and the line Hell 
in 8236 A. Lines of Nl and Nail between 8185 and 8242 A jRudv et al.l l2003h can also be 
found. 

Massive stars, suc h as some Wolf-Rayet, can exhibit emission lines of Pa9 (9229 A) and 
in Hell+Civ (8236 A) JCrowther eraDl2002h . 



The equivalent widths of the lines in the range covered by the narrow-band filters is 
of the order of tens of Angstroms. Although the relative frequency of appearance of these 
types of stars is very low, they can be also selected by their color excess. 



3.1. Classification schemes 



Further analysis is needed in other to classify the different objects selected by their 
narrow-band flux excess. 

To separate the stars, we can appl y morphological rn e thods to the selected candidates. 
A example of that approach appears in lMcCracken et al.l (120031 ). An alternative, and com- 
monly used, approximation is to ma ke use of the artificial n eural network of the source- 
characterization program SExtractor (IBertin fc Arnoutslll996l ). Nevertheless, these classifi- 
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cation methods depend on the SNR and seeing of the images. It can be difficult to classify 
candidates morphologically when the overall SNR is low or the seeing is not optimal. Mul- 
ticolor photometry can be also used to distinguish stars using a color criterion (we show a 
simple example in Sect. 14.3. ip 



Methods to classify extragalactic sources range from statistical corrections (IJones fc Bland-Hawthorn 



20011) to color-color criteria and photometric re dshift determination. For example , iPascual et al.l 
(I2OOII ) recomputed the statistical correction of I Jones fc Bland-HawthornI (I2OOII ) to find that 
the estimated contamination from galaxies at z > 0.24 was negligible. 



(IFuiita et al.l 120031) uses thre e photometric bands {B, Rc and Ic) and the GISSEL96 
models (IBruzual &: CharlotI |2003| ) to develop a color criterion that allow the authors to 
select Ha emitters at z=0.24:. Fo llowing a simila r procedure, but using spectra of emission- 
line galaxies rather than models, iLy et al.l (120061 ) uses five photometric bands {B, Rq, V, i' 
and z'). 

Lya surveys use dropout technique in filters bluer than i in order to assure that the 



Lyma n break (A = 912A) lies bluewards the narrow-band filter. For example. IXaniguchi et al. 



(I2OO5I ) request that the candidates are not detected in B, V and R, in addition to other 
criteria color criteria. 



Finally, we mention a more complex clas sification scheme . It is presented in IWolf et al. 



(I2OOII ) and it is used in the CADIS survey (IHippelein et al.l l2003l ) . Twelve medium pass 
band filters are used, in addition to broad-band B, R, I, J and K', to classify the objects 
in three broad categories: stars, quasars (including AGN) and galaxies. 



4. Measured flux and equivalent width of the emission line 

Fluxes of the objects through the filters are obtained integrating the spectrum of the 
object through the filter profile. The narrow-band filter technique allows, not only to select 
the objects with a possible emission line, but also to compute the value of the line fiux and 
the equivalent width of the line. 

In order to obtain the line and continuum fiux for the selected objects, several simplifi- 
cations, that are studied in this section, need to be made. First, the line profile is assumed 
infinitely thin when compared with the width of the narrow-band filter. Second, the posi- 
tioning of the emission within the transmittance profile of the used filters is very close to the 
center of the narrow-band filter. With simple assumption we can obtain a mean wavelength 
for the emission line. Finally, the continuum fiux can be modeled by an analytic function. 
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either a power-law or a polynomial. It is critical to have a good estimate of the continuum 
flux. 



4.1. The infinitely thin line approximation 

The FWHM of the lines of star-forming galaxies is related with the mass of the object 
and typically is less than ~ 10 A. To estimate the influence of the finite width of the emission 
line we have integrated Gaussian profiles through Gaussian and rectangular narrow-band 
filters and computed the recovered line flux. The FWHM increases with redshift a factor 
1 + z, so the infinitely thin line approximation can be valid, recovering more than 80% of 
the line flux, up to z ~ 4 for very narrow filters (A ~ 50A). For wider filters, (A ~ 150A), 
we recover more than 80% of the flux up to 2; ~ 10. These limit redshifts are lower for 
sources with broader lines, such as AGN. The ELGs can still be selected using the methods 
described in Sect. EJbut the fluxes of the lines computed with this section's equations will 
underestimate the total flux. 



4.2. Line flux and EW of a source with an infinitely thin line 

We assume that the SED of a source at redshift z with a emission line of rest wavelength 
Xr can be expressed as the sum of a continuum contribution and an infinitely thin emission 
line with flux /i. The emission line is represented analytically by a delta of Dirac centered 
in = Ar(l + z). 

In the following we assume that we have identified the emission line detected inside 
the narrow-band filter. The exact wavelength of the line is also unknown, but it is located 
in a very narrow range inside the narrow-band filter. We handle the determination of an 
approximate redshift for the line in Sect. 14.2.21 

With the assumed conditions, the resulting SED is: 

/(A) = /e(A) + /i 6{X - A,) (4) 

with /i constant and /c(A) an arbitrary function of the wavelength. Note that a negative 
value of /i produces an absorption line. 

We define the mean fiux through the filter as: 

//(A)T(A)AdA _ //(A)T(A)AdA 
^ /T(A)AdA /T(A)AdA ^ ' 
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The mean flux of tlie SED of Eq. H] is: 



/ = /c+('^V^)/i (6) 



The mean flux is a function of the redshift of the source z (or a function of A^, as both 
are related)E. The mean flux depends on the normalized transmittance of the filter at A2, 
the wavelength where the emission line is located and on the mean value of the continuum 
flux. We have not introduced yet any constraint to the shape of the continuum. 



4-2.1. Effective width of the filters 
Returning to Eq. [6l we introduce the effective width of the filter A': 

A'(AJ . (^) (^) (7) 



This definition is slightly different of the width obtained integrating the normalized 
transmittance. Both A and A' are independent of the maximum transmittance of the filter. 
A' depends on the redshift of the source. 

The value of A' is close to A in the center of the filter, in the region where the relative 
transmittance is close to unity. In the wings of the filter, the effective width increases and 
tends to infinity when the transmittance tends to zero. The physical meaning is that a 
emission line does not contribute to the integrated flux if it is positioned on a wavelength 
where the transmittance is near zero. The term in Eq [6] containing the line flux is effectively 
zero when the effective width tends to infinity. 



The result of I Walled (jl990l ) is based on square filters. We can recover this result making 
the transmittance equal to unity and making Aq = A^. In that case, the effective width of 
the filter reduces to: 



In general, the flux equation [6] can be rewritten finally as: 



^We will use both variables interchangeably, as appropriated 
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4-2.2. Estimation of mean redshift 

We focus on the position of the emission hne within the narrow-band filter. The defini- 
tion of effective width depends on the unknown redshift of the source. As the candidates are 
selected by their color excess in a narrow-band filter, the object with emission lines will be 
distributed in a small range of wavelengths, given mainly by the transmittance of the narrow 
band filter. If dn is the number of galaxies with a emission line, whose rest wavelength is Ar, 
the number of selected objects is P{\z,mB,m-!<s) dn. If we suppose that the probability of 
selecting a given object can be factorized in a term depending on and in terms depending 
on the magnitudes, the mean value of an arbitrary function h{Xz) can be written as: 

- _ J h{\,) p{\,) dn ^^^^ 
J p(A^) dn 

We can assume also that the probability p{Xz) is dominated by the shape of the narrow- 
band filter and write: 

p{Xz) oc An ' (11) 

Assuming a constant comoving density of galaxies inside the narrow-band filter, the number 
of galaxies is proport ional to th e comoving volume V (for the explicit expression of the 



comoving volume, see lHogglll999l . and references therein). The final equation is: 



_ J hjXz) A'^\Xz) dn _ Jh{Xz)XzMXz)dV{z) 

jA'^-\Xz)dn JXzMXz)dV{z) ^ ^ 

The mean value of A^ can be directly computed from last equation. The mean values com- 
puted with this method depend slightly on the cosmology parametrization selected (except 
on the value of the Hubble constant) and on the emission line chosen. 



4.3. Line flux and EW with different filter layouts 

The narrow-band filter layout used determines the approach used to approximate the 
continuum near the emission line. In a general case, we have n filters and n equations like 
Eq. [6l one for each filter. This provides us with, at most, n — 1 parameters to parametrize 
the continuum and one parameter to characterize the line, (/i). 

A power-law continuum (/^ oc A'^) can be fitted with at least three bands. It is even 
possible to assuna e a given power-law index (for example, the /3 = — 1 for Lya emitters in 
Ando et al. ( 20od )). Another possible functional form for the continuum is a polynomial of 



degree n — 2. 
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Most surveys use one narrow-band filter on the region of interest an d one or more 



broad-band filters (e.g. iPascual et al.ll200ll : iHu et al.ll2004l : lAjiki et al.ll2006al ) . Other surveys 



use more several co ntiguous narrow-band filters and a medium-band fi lter as a veto filter 
( IWestra et al.ll2006l ) or more complicated layouts (IHayes fc Ostlinll2006l ). 

In the following sections we study two different cases of the polynomial functional form. 
With two broad-band filters, a linear dependency of the continuum fiux with wavelength 
can be assumed. In the case of having one broad-band filter, the continuum has a constant 
value. 



4-3.1. Three filters and a linear continuum 

We develop here a general solution for the case of three filters covering the line. Solutions 
for other cases with more filters can be easily obtained. 

The continuum is approximated by a linear term with slope 7. For simplicity, we select 
the parameters so that the linear term vanish at A^: 

/c(A) = /c(A,)+7(A-A,) (13) 
Using Eq. [5], the effective flux is: 



/=/.+,(A„ + ^_AJ + _^ (14) 



being fi the second order central moment of the transmittance of the fllter (Eq. [3]). We 
rename the factor of 7 as s: 

s = Ao + ^ - A, (15) 

We also substitute A^ by A^ in s and in A'. Equation is written here for the three fllters 
covering the line (denoted A, B and N). 

/a = /c + 7SA + -^/i (16) 

^A 

/b = fc + lSB + ^fl (17) 
^B 

/n = /c + 7SN + -^/i (18) 

This is a system of three linear equations with three unknowns (/c, 7 and /i) that can be 
resolved easily. 
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4-3.2. Two filters and a constant continuum 

In this case, we denote the fihers by B and N and suppose that N is a narrow-band 
filter. 

The fluxes inside the fllters are: 

/b = /c + -^/l /n = /c + -T^/l (19) 

^B 

We account for the different effective width of the filters introducing the new parameter 
e = A^/Ag. As N is a narrow-band filter, e -C 1. The parameters of the emission line are 
as follows: 

/i = A^(/n-/b)-^ (20) 

1 — e 

/c = /b(^^) (21) 
W = A^(Q-1)(^^-^^ (22) 

being Q = /n//b- The color of a source (with the color normahzation of Sect. [2]) is directly 
rriB — rri^ = 2.5 log Q. 

Equation [22] shows that the EW does not change linearly with Q. There exists an upper 
limit to the color that a source can have due to the presence of an emission line: 



Q^ax = „lim QiW) = \im = -rf (23) 



It is worth noting that observing an emission-line object with a narrow-band filter 
increases the contrast between the emission line and the continuum, as the object is brighter 
in the line than in the continuum. Let us assume that we have two objects with the same 
magnitude in the broad-band, one with line of equivalent width W (named L) and the other 
without emission (denoted C). The ratio of the fiuxes of the two objects measured in th e 



narrow-band filter is the contrast parameter, defined originally by [Thompson et al.l (119951 ). 
Using Eq. [HI we obtain: 

/(L) 1 + 1^ 



/(C) 1 + 1 



1£L 
W_ 

t 

B 



(24) 



In the original contrast parameter of [Thompson et al.l (Il995l ). the width of the filter is char- 



acterized by its FWHM. For our Eq. [24] we, more accurately, use the effective width as 
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defined in Eq. [71 The numeric expression of the contrast parameter is equal to the equation 
obtained by inverting Eq. [221 

There are two hmit cases for these equations. When the width of the broad-band is 
much greater than the narrow-band width, we can consider that the hne do not contribute 
to the broad-band filter. The parameter e tends to zero and the equations can be written as: 



/i = A;,(/n-/b) (25) 
fc = h (26) 

W = A'Jk-l)=A'^(Q-l) (27) 



B 



The second limit case is when one filter do not cover the emission line. Thus, the line 
do not contribute to the flux in that band and the final equations are again Eqs. [251 - [271 



4.4. Multiple lines inside the filter 

In most cases of astrophysical interest, more than one line lies inside the narrow-band 
filter. For example. Ha plus [Nil] AA6548,6584, and [Olll] AA4959,5007 plus B.f3. The flux 
recovered from the narrow-band filter comes from all the lines, although each line experiences 
a different filter transmittance. 

We write Eq. [H with three different lines (the case of Ha and [Olll] AA4959,5007) and 
denote their fluxes and wavelength with super index 1, 2 and 3: 

/(A) = /e + f,' 6{X - Xl) + 5{\ - A^) + /3 5{\ - \l) (28) 
and Eq. [9l is substituted by: 

^^^'^ A^)^' + A^^' + MAf) ^^^^ 

The wavelengths of the three lines can be obtained from the redshift (or from the wavelength) 
of one of them (for example, A^ and A^ as a function of A^). We suppose that line 1 is the 
line we are interested on. 

If the emission lines fluxes have known ratios, we can simplify the equation: 



' ^'+T7^f) (30) 



) A'(A^)/ji A'(A3)/i 
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The effective widtli of tlie filter is substituted by a combined widtli (A ) tliat depends 

on tlie particular set of lines studied and their ratios, that we denote r2 = -k and = -k- 

Ji Ji 

The combined width is: 

A" (A!, r.. r,) . + + ^.3) ' (31) 

The equation developed previously to obtain fluxes and EW are valid to obtain directly the 
flux of the main line (line with index 1) using A" instead of A'. The mean wavelength is 
altered also. It is still computed using Eq. [T21 but substituting A' by A". 

The ratio of the combined width with the width computed in the single line case repre- 
sents the fraction of the total flux that comes from the main line. We denote this ratio by 

-HK.r.,r,) = ^^§f^ (32) 

The fraction (f) depends on the shape of the filter and the different line ratios of the lines. 
We can obtain a limit case of when the filters are wide (this is the case for broad-band 
filters). With that assumption, A'(A^) ^ A'(A^) ~ A'(A^) and Eq. [32 can be written as: 

0z(r2,r3) = (l + ra + rg)-^ (33) 

In the following subsections we apply Eq. [31] - [33] for the particular cases of Ha and 
[Oiii]. 



4.4.1. Ha and [Nu] 



Both lines of the doublet [Nil] are separated about 20 A from Ha in the rest frame. In 
this case, Aj is the redshifted wavelength of Ha, A^ is the redshift wavelength of [Nll]A6548 
and A^ is the redshifted wavelength of [Nll]A6584. 

The global ratio I([Nll])/I(Ha) is related with m etallicity (Ivan Zee et al. Ill998| : [Melbourne fc Salzerl 



2002 



1987 



Kewlev fc Dopi 



Swinbank et al. 



a 



2002l : [Pettini fc Pagelll2004l l and AGN activity (IVeilleux fc Osterbrock 



20041 ). showing a great range of variation between star- forming galax- 



ies. 



Historically, Ha surveys that could not separate [Nil] from Ha have used mean values 
to correct for the contribution of [Nil] . A commonly assu naed mean value o f the flux ratio 
I([Nll]A 6584)/I(Ha) is 0.32 . This is the value reported bv iKennicuttI fll992h and the UCM 
survey (iGallego et al.[ [19971 ) (I James et al.[[2005[ . see also). 
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This particular mean value is wi dely used in studies of the SFR based on Ha, where the 
lines of [Nil] can not be resolved (e.g. 



2OOOI : IPascual et al.ll200ll : iFuiita et al 



Tresse &: Maddox 



2003 



1998 



Umeda et al. 



20M). 



Yan et al. 1999: Iwamuro et al. 



The sar ne line ratio is obtained from t he Sloan Digital Sky Survey Data Release ^ 
(SDSS DR4, kdelman-McCarthv et aPbooeh if we apply a lower cutoff in EW of 20A (the 
detection limit of the UCM survey), typical of classical emission-line galaxy surveys. 

In Fig. [T]we show the ratio I([Nll])/I(Ha) represented as a function of the equivalent 
width of Ha for the galaxies of the UCM survey and the SDSS. The equivalent width rep- 
resented is observed, not rest frame. As the redshift of the main galaxy sample of SDSS is 
at most z ~0.1, the uncertainty produced by not correcting is about 4%. The effect in the 
UCM is even less, as the galaxies in the sample are local. 

The mean value assumed is represented as a straight line. We can see the wide range 
of variation of I([Nll])/I(Ha), from near 0.02 to 4. There is a tendency of decreasing 
I([Nll])/I(Ha) with increasing EW, corresponding to a decrease of the metallicity. 



^The intensity of the nitrogen lines is related: I([Nll] A6584)/I([Nn]A6548) = 3 (jOsterbrock 

19891 ). 



We can write Eq. [3T] as a function of r=I([Nll]A6584)/I(Ha) only. In this case, the 
combined effective width of the filter is: 



A"(A» 



1 



+ 



3A'(A2) A'(A 



-1 



(34) 



We obtain the mean redshift of the sources using Eq. [T2] and inserting there the 
A"(A^,r). We can obtain also the limit parameter 0; = (1 + 4/3r)~^. With the mean 
value assumed of r = 0.32, 0^ = 0.71. This means that, in the ideal case of a square filter 
covering the three lines and the given line ratios, 71% of the flux recovered in Eq. [20] comes 
from Ha and the rest comes from [Nil] . 

A narrow-band filter with the same value than the broad-band filter is wide enough 
to contain simultaneously the two nitrogen lines and Ha. A greater value of implies 
that the [Nil] lines only enter partially inside the narrow line filter. Using the standard 
[Nil] AA6548, 6584 correction with very narrow filters produce a systematic under estimation 
of the line flux. 

We can further refine the [Nil] AA6548, 6584 correction. From Fig [1] we can obtain a 



http://www.mpa-garching.mpg.de/SDSS/DR4/ 
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relation between the equivalent width H^(Ha) and r=I([Nll]AA6548,6584)/I(Ha;). Using 
that relation we can iterate through Eq. [22] and obtain simultaneously r and EW. With this 
procedure, every object has a different value of r, which implies different values of A" and </> 
for every object. 



4.4.2. H(3 and /Olll] 

In this case, three lines enter the filter, the doublet [Olll] and H/?. H/3 is separated 146A 
from [Olll] A5007 in the rest frame. We can write Eq. [3T]with as the redshifted wavelength 
of H/9, A^ is the redshift wavelength of [Olll]A4959 and A^ is the redshifted wavelength of 
[Olll]A5007. 

In Fig [2] we show the ratio [Olll] A5007/H/3 again for galaxies of the SDSS DR4 and for 
galaxies of the UCM sample. We have only represented galaxies with the lines [Olll] AA4959,5007, 
H/? and a EW(Ha) > 20A, in order to mimic the limit values of classical ELG surveys. The 
mean value obtained for I(H/5)/I([0lll] A5007) in the SDSS is r = 1.05 (the result for the 
UCM is compatible, r = 1.1). 

We write Eq. [HI] for this particular function of r = I(H/3)/I([0lll] A5007) and 

knowing that I([Olll] A5007)/I([Olll]A4959) = 3 JOsterbrocklll989l ): 



1 1 1 



The limit parameter is in this case 0/ = (4/3 + r) ^. With the mean value assumed for 
r, (pi = 0.42, meaning that only 42% of the flux would come from [Olll]A5007. 

In this case, is very close to unity for the narrow-band filters. This means that 
[Olll]A5007 will enter alone in the filter. In the case of broad-band, the ratio is close to the 
limit value, meaning that the three lines enter the broad-band filter. As a summary, we can 
consider that [Olll]A5007 enters alone in the narrow-band filter, but we have to include the 
three lines when considering the broad-band filter. 



4.5. Line luminosity of the objects 

The line luminosity of each object can be calculated from the measured line flux and 
the redshift of the object (obtained from A^). The line flux corrected from obscuration by 
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^ -1.5 - 



-2 ' ' 1 1 1 ^LJ 

1 1.5 2 2.5 

log EW(Ha) [A] 

Fig. 1. — Ratio I([Nll]A6584)/I(Ha) represented as a function of the equivalent width of Ha 
for the galaxies of the UCM survey (white circles) and the Sloan Digital Sky Survey (the 
number of galaxies is represented in gray scale). The mean value assumed in this work is 
represented as a straight hne. 




-1.5 I 1 1 1 

0.5 1 1.5 2 2.5 

log EW(HP) [A] 

Fig. 2. — Ratio I(H/9)/I([0lll] A5007) represented as a function of the equivalent width of 
H/5 for the galaxies of the UCM survey (white circles) and the Sloan Digital Sky Survey (the 
number of galaxies is represented in gray scale). The mean value assumed in this work is 
represented as a straight line. 
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dust (/°) can be obtained from the measured line flux (/i): 

log/° = log/i + 0.4(AA„ + aAj (36) 

being the internal obscuration at the rest frame wavelength A,- and ax^ the galactic 
obscuration. The line luminosity is then: 

L = 47r dliz) /o (37) 

being di^{z) the luminosity distance. 

The common obscuration corrections applied in narrow-band searches are discussed in 
the following subsection. 



4-5.1. The obscuration correction 



Obscuration by dust affects seriously any measurement of galaxy luminosities at the 
optical and UV wavelengths. Correcting for this effect is a difficult task as measurements of 
parameters sensitive to obscuration are not always easy to obtain. 

In the case of obscuration corrections for emis sion lines, it i s coni mon to assume an 
average value of 1™ in Ha and the extinction curve of lCardelli et al.l (119891 ) to apply the value 
to other emission lines. This value of Afja has been fo und as the average of the obscuration 
of several samples o f local Ha emitting galaxies (e.g. lKennicuttlll983l : iGallego et al.lll995 



Hopkins et al 



20031 ) . This procedure has the advantage that is very simple to apply to a 



luminosity function or luminosity density, since it is a scaling factor. 

It is worth to note here that the average o bscuration applied in UV studies is about 
At/ = 1"" hresse et al.lE996l : ISuUivan et allboool . e.g.). In the GALEX far-ultraviolet (FUV), 
Buat et al.l (120051 ) found a mean obscuration of Apuv = 1-1™ for galaxies selected in the 



GALEX near-ultraviolet band. This value is roughly Ay 0.4"^. Galaxies in the GALEX UV 



Imaging Atlas of Nearby Galaxies (IGil de Paz et al.l |2004| . |2006| ) have a mean obscuration, 
considering only spirals and irregulars, of Apuy = 1.8™", corresponding to Ay ~ 0.7™ (Gil 
de Paz, private communication) as this sample can be considered as selected in the optical. 

The average obscura tion of Aho = 1™ corresponds to a value of about Ay ^ 1.22™ as 
noted by iHopkind (120041 ). This disagreement in the obscuration is a selection effect of the 
different (FUV, UV and Ha) samples. Galaxies selected at redder wavelengths can have 
greater internal obscurations than those selected at the UV, where heavily absorbed systems 
will not be selected. 
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The value of A at the wavelength of the emission line is obtained applying a given 
extinction curve. We are assuming that the continuum and the emission line are affected by 
the same obscuration, so the EW does not ne ed obscuration correc tion. This last assumption 
is not true if we use the obscuration recipe of ICalzetti et al.l (Il996l ). that is applied for strong 
starbursts. However, normal galaxies and moderate starbursts can be assumed as showing 
no obscuration correction in EW. 



5. Magnitude-color selection criterion 

In this Section, we study the underlying statistics of the magnitude-color diagram used 
to select ELGs in the most simple scenario, a broad-band filter and a narrow-band filter. 
A sample selected in broad-narrow color is directly selected in EW (Eq. [22]). Thus, ideally, 
a sample can be created selecting the objects with color in excess of a minimum color 
(me — fnj<f)^i^ that translates into a minimum equivalent width EWmin. 

The dispersion in the magnitude-color diagram is dominated by the uncertainties in 
the magnitudes in the faint region and by the different colors of the objects in the bright 
region. Furthermore, the mean value of the color, /i[mB — mj<{], depends on mN, because 
different types of objects can dominate the color at different narrow-band magnitude ranges. 
The method of selection has to account for the different sources of dispersion of the color. A 
solution is to compute the standard deviation of the color dispersion and consider candidates, 
with a level of signification n^, the objects with color: 

me — > max{(mB — 'mN)min, At["^B — "^n] + naCrlrriB — TTiN]} (38) 

The standard deviation of the color distribution, as a function of the narrow-band 
magnitude, can be calculated from the distribution of objects. In the following subsection, 
we infer an analytic expression that can be useful when there are not enough points in the 
data to calculate directly the dispersion. 



5.1. Standard deviation of the color distribution 

Let us assume a population of objects with the same SED, given by /(A), so all the 
objects have the same color. 

In general, the source flux in analog to digital units (ADU) collected in a CCD detector 
from an object are proportional to the integrated flux of the source through the pass band 
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and the exposure time t: 

N (xtj /(A)T(A)AdA (39) 

Equation [39] can be written in terms of Aq and A as follow^: 

AToct/TpeakAAo (40) 
where Tpeak is the maximum value of the transmittance. 

The ratio Q = /n//b is proportional to the ratio of ADUs in the same bands TZ: 

^^^.am(^)(^M(hi).,aM (41) 



In general, Q is a function with a strong dependency on the wavelength because of night sky 
emission lines. In the cases covered in the present work, the filters are selected to exclude 
strong emission features in the sky spectrum but that is not always possible, such as when 
using narrow band IR filters designed to cover emission lines of astrophysical interest at zero 
redshift like Br7 or molecular hydrogen. 

We include the parameters depending of the transmittance of the filters in: 

(42) 




Eq. HI] applies independently of the flux source. We denote the quantities referred to 
the sky background with a index S and the quantities referred to a particular object by O. 
In particular: 



The standard deviation of Qo can be obtained from the standard deviation of the ratio 
of the source flux in ADUs: TZq as: 



Qo J [no J { NO J ^[ NO 



(44) 



Equation Ho] can be also written in terms of Acff ( Fiorucci fc Munari 20031 Eq. 24) as follows: 



N CCtf Tpeak A Acff 

being f-imiWdx ^nd A ff - f^f^rwdx 
Demg / - ^ ^^^^ ana Aeff - ^ ^^^^ ^^^^ 

We prefer using Aq to avoid the flux dependence of Aofr, that would complicate the equations of Section 21 
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We do not include a covariance term. The covariance between the sources flux in ADUs 
should be negligible because we are measuring different parts of the spectrum with very 
different filters. 

The terms in the right-hand side of the equation are the inve rses of the signal-no ise ratio 



(SNR) We adopt here the signal to noise ratio as computed by iHowell et al.l (119961 ) with a 
slightly changed notation: 

SNR ^ ^ = = (45) 

Nq are the ADUs of the object in a given band, the best estimation of the sky background 
per pixel, A the area in pixels of the aperture used to measure photometry, R is the readout 
noise in ADUs and G the inverse of the gain g. The term involving the readout noise can be 
neglected when the exposures are long enough to be background limited, so the dominating 
contribution to the noise are the Poisson terms. Note that the background limit can be 
difficult to reach if the narrow-band filter is particularly narrow, as the airglow windows 
considered in the present work are regions of low sky background. 

Equation HS] depends on the ADUs of the object and on the sky background measured 
in both bands. We use Eq. HTJto translate the flux in the broad band into flux in the narrow 
band, assuming a given color for the object (Qo) and for the sky background (Qs)- Finally, 
we put the equation as a function of the ADUs per unit time and pixel of the sky (ra^) and 
of the ADUs per unit time of the object {N§), so we can easily operate in magnitudes: 

^[QoW _ + 



Qo J Vn tsj N§ 



+ ^G^flM;^ + t^S6;^)l^) (46) 



We introduce a new parameter q that indicates if we are in the noise regime dominated 
by the background or in the regime dominated by the Poisson noise produced exclusively by 
the source flux. The parameter q is the ratio between the flux in ADUs (or ADUs per unit 
time) coming from the sky background and that coming from the object, inside the given 
aperture A and measured with a certain filter: 



(47) 
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where n is the surface brightness of the sky, m the magnitude of the object and p the plate 
scale, q tends to zero when the background level is low compared with the brightness of the 
object and to infinity when the background dominates over the object. 

The previous equation is written in terms of gN, but can be rewritten easily in terms of 
ge, as gN = Qb Qs/Qo- 

Equation US] shows (t[Qo] as a function of N^, and through it, as a function of the 
narrow-band magnitude. Of the three terms of the equation, the last corresponds to the 
readout noise and is negligible except for very short exposure. The second term comes from 
the contribution of the sky background and dominates for faint objects. The first term is 
produced by the ADUs of the object and only is important (when compared with the sky 
contribution) for bright objects. 

The presence of a logarithm in the transformation between flux and magnitude produces 
an asymmetric distribution of points around the mean color, mainly in the region of faint 
objects. We obtain this behavior by defining two curves 



with the level over the mean. 

As can be seen from Eq. HBl for a given narrow-band magnitude, cr[Qo] depends on 
the properties of the filters encoded in b, on the exposure times in the two bands {t-^ and 
te), on the color of the sky background (Qs) and on the relative brightness of the sky when 
compared with the brightness of the objects (g). In Fig. [3] we show how the selection curve 
changes (with a fixed value of n„=3) with the exposure time and sky background in the 
two bands for the sample filter set CAFOS8200 (the different filter properties are detailed in 
Sect. ED The increase of exposure time or the reduction of the sky background (effectively 
increasing the sky-background magnitude) reduces the slope of the selection curve. 



Isolines represent curves joining points with an equal given property. Lines with equal 
EW correspond with lines of equal color in the magnitude-color diagram, but lines with equal 
continuum or line flux are not directly related with neither the color nor the magnitude of 
any of the bands. 



(me - mw)-; 
(mB - mN)+; 



2.51og(Q-n<,a[Q]) 
2.51og(Q + r2,a[Q]) 



(48) 
(49) 



5.2. Isolines in the magnitude-color space 
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Fig. 3. — Selection curves with several different exposure times (top) and sky background 
(bottom). The broad-band exposure times are: 200, 500, 800 and 1100 s, whereas the 
exposure time in the narrow-band filter is 800 s. The sky background broad-band magnitudes 
are 17^", 17.5™", IS™ and 18.5"*. We represent also the 5-a broad-band limit magnitude (lines 
crossing on right top) corresponding to each exposure time and sky background magnitude. 
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5.2.1. Isolines of line flux 



We can compute the isolines in magnitude-color space starting with Eq. [91 If we assume 
a constant line flux f\, the magnitudes of the objects as a function of f\ and EW only are: 



TUB 



niB — TTiN 




(50) 
(51) 

(52) 



The luminosity of the emission line is proportional to the measured flux (see Sect. 14.51 for 
details), so the isolines of line flux are also isolines of line luminosity. 

We can write a parametric expression for the isolines of line flux. We use w = W/ A'j^ as 
an adimensional parameter. Note that the equations are also valid for absorption lines (/i< 
and W < 0). To simplify the notation, we rename the color tub — ttin as the independent 
variable y and the narrow-band magnitude as X . 

The parametric equations of the line flux isolines are: 



X = a;oo-2.51og 
2.51og' 



1 

1 + - 

w 



1 + ew 



(53) 
(54) 



where Xoo 
--C - 



X 



is the narrow-band magnitude of an object dominated by the emission line: 
2.5 log -jn- ■ Equations [53] and [Ml are valid in the ranges of emission lines, cor- 
responding to w G (0, oo), and in the domain of absorptions lines, where the adimensional 
parameter w G (—1,0). 

The asymptotic behavior of the curve is as follows. For an object without line w 0; 
the color ?/ — 0. When w increases asymptotically the object is dominated by the emission 
line. Then , x — > x^o (the magnitude of an object dominated by the emission line) and the 
color y — > — 2.51oge. When w — > —1, the object is dominated by the absorption line; the 
absorption line grows stronger and the narrow-band magnitude tends to positive infinity; 
the object is progressively dimmer as the absorption line increases its equivalent width. The 
broad-narrow color tends to negative infinity. 
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5.2.2. Isolines of continuum flux 



The isolines of continuum flux can also be obtained 

rriB 



C- 2.5 log 
C- 2.5 log 



W 



(55) 
(56) 



The equation of the color me — is equal to Eq. [521 because the color depends only on 
the equivalent width. A similar approach to obtain the equation of the isolines of continuum 
flux can be used. The adimensional parameter is again w = W/A'-j^ and the parametric 
equations are: 



X = xo — 2.5 log (1 + u;) 

, , l + w 
y = 2.5 log 



1 + ew ^ 

being xq = C — 2.5 log (/c) the magnitude on an object without line and w G (—1, oo). The 
asymptotic limits correspond to a object without line {w —>■ 0, x xq and y 0), an 
object dominated by the emission line (w — > oo, x — » — oo and y — > —2.5 loge) and an object 
dominated by an absorption line {w —>■ —1, a; ^ oo and y — oo) 

In Fig. m we show a magnitude-color selection diagram including example isolines for 
the WFC9200 case (see Sect. [6] for details). The isolines of continuum flux and line flux are 
represented by the dashed and solid curves respectively. We can see clearly the asymptotic 
behavior of both types of isolines. Although the limited numeric precision of the plot do not 
show it, the curves of line flux continue up to the maximum color. 



(57) 
(58) 



5.3. Redshift range of detectability and covered volume 

We have assigned an approximate redshift z to the source, but we do not know which 
is the redshift range Zi z^ of the objects that could be selected and included in the sample. 
The answer to this question is the key for the incompleteness correction. If we intend to 
compute a luminosity function, we need to compute the comoving volume per object. 

We consider that the object is at redshift z and it has a line flux f\ and continuum 
flux /c. The flux recovered in a given filter at redshift z (whereby the emission line lies at 
wavelength ) is: 
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Fig. 4. — Isolines of line flux (solid curve) and continuum flux (dashed curve) for emission 
lines (broad-narrow color above zero) and absorption lines (bellow zero) in images obtained 
with the WFC9200 system (see Sect. [6] for details). The flux lines correspond to 10~^^ and 



10 erg s ^ . The continuum lines are 10 



16 



10-1^ and 10-is erg s~i A~i. 
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with (11 the luminosity distance. For 2; = z we recover Eq. [91 Substituting the appropriate 
filter widths A', we can obtain the magnitude and color of the object as a function of redshift. 
Only at a range of redshifts Zi, Z2 the object fulfills the color criterion given in Sect. O 

In Fig. [5] we show the redshift traces given by Eq. [59]for two objects in a magnitude-color 
diagram corresponding to the filter set CAFOS8200 (see Sect.[6]for details). The position of 
the objects in the diagram at redshift z are represented as crosses on top of the traces. The 
intersection between the traces and the selection curve (drawn with a dashed line) is labeled 
as zi and ^2- The object would be selected only in the range zi < z < Z2. 

The values of zi and Z2 depend through Eq. [591 on the transmittance of the filters. 
As the transmittance of the broad-band filter is mainly constant in the wavelength range 
of interest, the narrow-band filter dominates. Thus, the redshift range Zi, Z2 is small and 
dominated by the shape of the transmittance of the narrow-band filter. The redshift range 
depends also on the selection curve. Objects selected at lower signification level are closer 
to the selection limit and the redshift range zi, Z2 is smaller. 

The comoving volume where the object could be selected is computed from the redshifts 
Zi, Z2. This volume is one of the ingredients needed to compute the luminosity function. With 
Eq. [59land the selection curve, we can compute a comoving volume for each object. It can be 
seen clearly that, with the same EW, faint objects have smaller comoving volumes. These 
objects have, for a fixed EW, low line flux (and low line luminosity). Thus, assuming a 
global comoving volume for all the detected objects tends to underestimate the value of the 
luminosity function for the faint objects and objects selected near the selection curve. 



6. Example filters 



6.1. Sets of filters and instruments 



In order to provide observational tests for our analysis, we select three sets of filters in 
two different telescopes. The first instrument is CAFOS on the 2.2 m telescope at Calar Alto 
Observatory (Almeria, Spain). The filters used were 816/16 (narrow band) and 850 /150c 
(broad band, Johnson I). This is the instrumental setup used in the [Pascual et al.l (1200 ll ) 



survey. The second instrument is the Wide Field Camera (WFC, Ives et al. 19961 ). located 



at the Isaac Newton Telescope in the Observatorio del Roque de los Muchachos (island of La 
Palma, Spain). Two different narrow-band filters have been used, #209 and #208, together 
with the broad band filter #194 (RGO /). 



The transmittance of the different filters is shown in Fig. [61 The fiducial transmittances 
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17 17.5 18 18.5 19 19.5 20 20.5 21 



Fig. 5. — Redshift traces (solid lines) of two objects marked as crosses in a magnitude-color 
diagram for CAFOS8200 (see Sect. [6] for details). The emission line in this case is Ha. The 
selection limit is the dashed line. The intersection between the traces and the selection curve 
delimits the redshift range (from zi to Z2 in the plot) where the object could still be selected. 
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are displayed together with the QE of the detector and the total transmittance. For a working 
temperature of 0°C, the blue shift produced by the temperature and the converging beam for 
CAFOS8200 is lOA, about 5% of the width of the filter (S. Pedraz, private comr nunication) 



The blue shift is 17 A for the filter WFC8200 and 19 A for the filter WFC9200 flKing et al. 



19941 ) ■ These two shifts are about ~35% of the width of the filters. The converging beam 
can produce also an effect in the absolute transmittances and widths, but it turned out to 
be negligible for our examples. 

Several quantities useful to describe filters are summarized in Table [U Columns 1 and 
2 are the name of the filter set and the name of the filter (either broad- or narrow-band). In 
the next three columns we show Aq (column 3), \fj^ (column 4) and A (column 5) computed 
directly from the transmittance. Column 6 is the estimated mean value of computed with 
Eq. [121 concordance cosmology and the Ha line, not including the corrections for the [Nil] 
lines (see bellow). With we compute the effective width of the filter in column 7. The 
difference between commonly used cosmologies and different emission lines is less than 5% 
in all cases. Note that the filter #194, when used together with #208 has an effective width 
roughly three times its width. The wavelength where the emission line is located has low 
transmittance. In this case, using A instead of A' would underestimate the continuum flux. 

In Table |2l we summarize the mean values of the combined width and parameter 
(Eq [32] for Ha and [Olll] -computed quantities with concordance cosmology. Columns 1 
and 2 have the same meaning in Table [H Columns 3 to 5 correspond to Ha-computed 
quantities with r =I( [Nil] A6584)/I(Ha) =0.32 (Sect. WA^\ and columns 6 to 8 to [Olll] 
with r =I(H/?) /I([Olll] A5007)=1.05 (Sect J4.4.2]) . Columns 3 and 6 are the mean wavelengths 
computed according to Eq [121 using the value of A" . Columns 4 and 7 the value of the 
combined width at the computed mean wavelength value and columns 5 and 8 the parameter 
0, that denotes the the fraction of the total flux that comes from the main line. 

For Ha, the narrow-band filter in the CAFOS8200 filter set has the same value than 
the broad-band filters. This means that is wide enough to contain simultaneously the two 
nitrogen lines and Ha. The narrow-band filters of WFC8200 and WFC9200 have a greater 
0, meaning that the [Nil] lines only enter partially inside the narrow line filter. In this two 
last cases, 82% and 85% of the flux in the narrow-band filter come from Ha, instead of the 
nominal 71%. Using the standard [Nil] AA6548, 6584 correction with very narrow filters such 
as the narrow-band filters of WFC, produce a systematic under estimation of the line fiux 
of about 15%. 

In the case of [Olll], is very close to unity for the narrow-band filters. This means 
that [Olll]A5007 will enter alone in the filter. In the case of broad-band, the ratio is close 
to the limit value (0.42 with the assumed value of r), meaning that the three lines enter 



1 



0.8 - 




7000 7500 8000 8500 9000 9500 10000 



Wavelength [A] 

Fig. 6. — Transmittance curve for the different sets of filters. Dasli-dotted lines represent the 
quantum efficiency, dotted lines represent the fiducial transmittance of the filters. Finally, 
the solid lines represent the real transmittance corrected from the different contributions 
detailed in the text. On top, the filters from CAFOS. Broad-band is 850/150c and narrow- 
band 816/16. On bottom the ffiters from WFC. Broad-band is ^^194 and narrow-band are 
#209 around 8200A and #208 around 9200A. 
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the broad-band filter. As a summary, we can consider that [Olll]A5007 enters alone in the 
narrow-band filter, but we have to include the three lines in the broad-band filter. 



6.2. Colors of galaxies 

We use galaxy templates to calculate th e color evolution of d ifferent galaxies with in- 



creasing redshift. With the set of templates of iKinney et al.l (119961 ) . we have obtained colors 



for the three groups of sample filters. The templates have a resolution of SA. We compare 
three wide types of galaxies with similar colors in this particular system: ellipticals and 
bulges, early-type spirals are represented by a Sb and finally late-type spirals and starbursts 
are represented by the starburst template SBl. The Sa and Sb spectra contain the Ha line 
in emission. The Sc and the starburst spectra also contain the nebular lines [Oll]A3727, H/3 
and [Olll]AA4959,5007. 

We show the color uib — '?7^n of the different templates as a function of redshift up to 
z = 2 in Fig. [7] for the three filter sets. Only the effect of the redshift is considered. We 
are not including galactic evolution. As expected, galaxies with different emission lines show 
strong increases on their colors when the emission lines enter the region covered by the filter. 
In the filter set WFC9200, elliptical templates exhibit color excess at redshifts around 2; ~ 1 
and higher. This excess is not produced by emission lines. Instead, the very different central 
wavelengths of the two filters involved and the slope of the continuum near the 4000A break 
explains the broad-narrow color obtained. 

An additional color can help to distinguish different types of objects with narrow-band 
color excess. Figure M shows a / - NB versus I — Z color-color diagram for the WFC8200 
filter set. We represents simultaneously stars, the black-body and galaxy templates. With a 
/ - NB versus NB diagram we can not distinguish between the color excess produced by the 
emission line of a starburst galaxy and by the molecular bands of a late-type star. But these 
objects have different I — Z color. Late type stars are very red, I — Z > 0.4 and starbursts 
are blue, with I — Z < 0.4. Thus, we can classify objects using their I — Z color. 



6.3. Colors of stars 



Using the library of Pickles ( 1998 1^ we have estimated the colors of stars of different 
spectral types in our sample sets. The library contains spectral energy distributions of stars 



'''http: / / www.ifa.hawaii.edu/users/pickles / A JP /hilib.html 
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Table 1. Summary of properties of the filters used 

















System 


Filter 


Ao 




A 




A'(A^) 






(A) 


(A) 


(A) 


(A) 


(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


CAFOS8200 


816/16 (narrow) 


8139 


80 


173 


8147 


176 




850/150C (broad) 


8715 


643 


1665 


8147 


1783 


WFC8200 


#209 (narrow) 


8157 


17 


48 


8157 


48 




#194 (broad) 


8007 


477 


1333 


8157 


1940 


WFC9200 


#208 (narrow) 


9190 


17 


49 


9190 


49 




#194 (broad) 


8007 


477 


1333 


9190 


4398 



Table 2. Mean combined effective widths of the filters 
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Fig. 7. — Colors of galaxy SEDs for the sample filter sets as a function of redshift . From top to 



bottom CAFOS8200, WFC8200 and WFC9200. Hubble types templates from lKinnev et al. 



(Il996l ) are plotted. Ellipticals (and bulges) are represented by dashed lines, we use Sb 
template for early type spirals and SO (dash-dotted lines) and SBl starbursts template for 
starburst and late spirals (Sc) (solid line). The position of the different emission lines is 
labeled. 
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Fig. 8. — / - NB versus I — Z color-color diagram for the WFC8200 filter set. The key for 
the different types is represented in the box. Starbursts and late type stars exhibit color 
excess in /-NB but have very different colors in 1-Z. 
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with a large range in metal content. The spectral range covers up to 10000 A. 

In Fig. [HI the color — as a function of spectral type is plotted. As noted before, 
stars follow the trend the colors of are dominated by bands of TiO and VO. 

Comparing the colors obtained with CAFOS8200 and WFC8200 we can see that cold 
stars with molecular bands can be misclassified as ELGs in the case of WFC8200 and 
WFC9200. 



7. Optimal exposure times 

A critical question that arises when planning observations with narrow-band filters is 
how long we should expose in the narrow band filter in order to obtain images that can be 
compared with the corresponding broad band. Exposure simulators available in telescopes 
usually only cope with broad-band filters. 

To be selected in a magnitude color diagram with a signification ricr, an object with 
an emission line has to show color excess over the mean distribution of objects at a given 
magnitude m^. The ratio of fluxes Qi of an object with emission line has to fulfill: 

Qi> f^[Qo]+n^(r[Qo] (60) 

being Qo the ratio of fluxes of the objects without emission line and the same narrow-band 
magnitude that the object with the emission line. 

We can approximate /i[Qo] ~ Qo and then: 

±fa_A>^ (61) 

n, VSo y So ' 

If we include the SNR via Eq. HH and Eq. HS] we can obtain the limit condition for the 
signal-noise ratios: 

2 



/ = — =SNR-^ + SNR-^ (62) 

The SNR depends on the magnitude of the object, the sky background, the aperture where 
the photometry is measured, some parameters of detector (gain and readout noise) and in 
the exposure time. In our problem the only free parameters are the exposure times, as the 
others are fixed by the observing conditions and by the particular emission line object that 
we wish to detect. 



We intend to make the dependence in exposure time explicit. For exposures where the 
readout noise of the detector is negligible, the signal-noise ratio goes roughly as the square 
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Fig. 9. — Colors of stars of late spectral types in our three sample filters. Stars of late spectral 
types are also represented with triangles, squares and circles for CAFOS8200, WFC8200 and 
WFC9200. The colors of a black body on the same filters have been added for comparison and 
are represented by solid, dashed and dotted lines respectively. The bottom axis represents 
the effective temperature of the star obtained from the spectral type of the model (top axis) . 
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root of the exposure time when the other free parameters are fixed. For the broad-band, we 
comprise in the constant a all the terms not depending on time: 

SNRb = a V^b (63) 

The signal-noise ratio in both bands is related through the ratio of exposure times. If 
we start from Eq. HSl we can easily obtain (without any approximation): 



For long exposures in both bands, the terms including i?^ can be neglected and we 
obtain simply: 



Vsnrb; 



/5^ = 6QM^±^) (66) 



being /3 the following 

The parameter /? has limit cases depending on ge- If — ^ 0, the noise is dominated by 
the Poisson noise of the source flux: 

/3' = &Q (67) 

If ge ^ oo, the equation is dominated by the Poisson noise of the sky background. This is 
the case for faint objects in the images: 

^2 = 6^ (68) 



s 

We also make explicit the dependence on exposure time of the narrow-band combining 
Eq. [63] and Eq. [65] to obtain: 

SNRn = a /? (69) 

Equation [62] can be written: 

t-Q jj^ 



We try to detect 5-cr over the mean objects with ttin = 22 and Qi = 1.8 (this represents 
a color excess of 0.6 magnitudes ). Assuming a mean color of the objects without emission 
line of Qo = 1, we immediately compute p = 0.16. 
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io^ 10^ lO"^ 

tg [s] 

Fig. 10. — Exposure times in broad- (t^) and narrow-bands (t-^) that allow the detection 
of an object with a given color ratio Q\. The corresponding SNR is also represented. The 
different curves in solid line correspond to different values of the parameter p (see text for 
the definition). Prom left to right p = 0.16 (this value is used in examples in text), 0.14, 
0.12 and 0.10. The dashed line represents the condition of minimum total exposure times 
and the dash-dotted line the condition of equal SNR in both bands. 
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With these filters, the color of the sky background is Qs = 0.23 and b = 0.024. From the 
exposure time calculator of the INT WFC, SIGNAL^, we obtain that, for an object of 
= 22, with seeing of 1" and a brightness of the sky background of 19.60 mag/arcs^ (typical 
for a dark night), the SNR in the broad band is the square root of the exposure time times 
a = 0.7. 

For the narrow band we compute (3. In this case, qb, computed from the magnitude of 
the sky background and the magnitude of the objects without emission line, is ~ 300. This 
means that the noise is dominated by the contribution of the sky background and f3 = 0.32. 

In Figure [10] we represent the exposure times in both bands that fulfill the condition in 
Eq. [To] for different values of p, from 0.16 (the example we are considering) to 0.10 in steps 
of 0.02. For each value there is a set of different exposure times in both bands (given by the 
curves) that fulfill the conditions. 

Each curve is limited by two asymptotic limits, obtained from the condition of infinite 
SNR in one of the bands: 

FortN^oo 4^ (71) 



{pa] 



FortB^oo 4- (^2) 

In the current example, the limit exposures are tg™ = 80 s and tjjj™ = 780 s. As the 
curves lay over their asymptotes, this exposure times are lower limits. 

The problem of selecting the optimal exposure times in order to optimize selection of 
objects with the narrow-band filter still exists at this point. We have to impose one additional 
condition to extract a pair ts, from Fig [TO] It seems natural to obtain the point that 
minimizes the total exposure time te + i^N- This condition is equivalent to minimize the sum 
of inverse squared SNRs with the constraint of constant + ^n. 



The exposure times are: 



- ^ n ^ (73) 
p^ p 

' i«B (74) 



'http: / /www. ing.iac.es/ds/signal/ 
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In terms of the signal-noise ratio, this last equation can be written as: 

SNR2 = (75) 



P 



SNR^ = ^-^ = /?SNR2 (76) 



For our particular example, the exposure times are: = 330 s and = 1030 s. The 
corresponding signal to noise ratios are SNRb = 13 and SNRn = 7. 

From the previous equations we can extract an important conclusion. As typically 
(3 < 1, the exposure times in the narrow band are greater than the exposure times in the 
broad band, but the SNR in the narrow band is lower than the SNR in the broad band. This 
can be a problem in some extreme cases. The narrow-band image is used normally to detect 
the objects. Hence, it is necessary to have enough signal-noise ratio in the narrow-band 
image to be able to detect the objects in it. Note that these exposure times do not account 
for the overheads proper of observations (dithering, readout time of the detector, etc.). 

Other conditions can be applied to Eq. [70] and different exposure times are obtained. If 
we wish to assure that the SNR is high simultaneously in both images, we can impose the 
condition of equal SNR to the objects without emission line. Both terms on the right hand 
side of Eq. [70] are equal. The exposure times are: 

tB = ^ (77) 

= 7 ^ = ^ (78) 

and the common signal-noise ratio is: 

SNR| = SNR2 = A (79) 
P 

If we apply the condition of equal SNR to our example, we obtain tB = 160 s and 
= 1560 s, whereas the SNR in both bands is SNRb = SNRn = 9. The total time in this 
case is 1720 s, that can be compared with the total exposure time of 1360 s in the optimal 
case. We have represented in Fig. [10] the lines that correspond to minimum exposure time 
(dashed) and equal SNR (dash-dotted). 

Other possibility is to work in the asymptote of long broad-band exposure times. With 
the availability of multicolor imaging surveys, deep broad-band public images can be used 
instead of observing simultaneously in different bands, including the narrow-band. In this 



-42- 



case, the flux limits are completely determined by the depth of the narrow-band image (in 
Eq.m SNRi^ < SNRn'). 

As we have shown, the condition of minimum exposure times provides high enough 
SNRs to carry out object detection in the narrow-band image. In the SNR obtained is not 
enough, other conditions can be imposed. The presented exposure times are only lower 
limits, obtained from the limit condition of Eq. [62l Longer exposure times, if possible, 
provide better SNRs in both bands that allow determining the line flux and EW of the 
emission line with lower uncertainties. 

The goals of the emission-line survey should be considered when deciding exposure 
times. Searches for relatively rare objects should survey relatively more area at the cost of 
getting lower SNR on the narrow-band filters. On the other side, very deep narrow-band 
survey can select candidates that will have very difficult follow up spectroscopy. In this last 
case, very deep narrow-band imaging is appropriated, in order to recover line fluxes with 
very low uncertainty. 

8. Summary and Conclusions 

The main results of this work follow: 

1. We discuss the contamination of a narrow-band survey, either by galaxies or stars 
and show different classification schemes that can be used to separate their different 
contributions. In general, multicolor photometry can be used to classify stars and 
ELGS at different redshifts. Morphological criteria can also help in the classification. 

2. For different narrow-band filters, we compute the mean colors of stars and galaxies. 
With all the sample filters, the galaxies with emission lines exhibit a color excess 
and they can be selected as noted by previous authors. Stars of different types can 
be selected depending on the shapes of the filters involved. A particular analysis is 
needed in each case. In general, late type stars have molecular bands in the regions 
studied that can be easily confused with emission lines. 

3. The line fiux and equivalent width are obtained with different assumptions about the 
continuum and different filter layouts, line We also confirm that the assumption of an 
infinitely thin line is appropriate for most physical cases. 

4. We study the case of several lines entering simultaneously in the narrow-band filter. 
This problem arises for Ha plus [Nil] AA6548,6584 and [Olll] AA4959,5007 plus H/3. In 
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each case we compute the correction due to the contribution of the other hnes. For 
very narrow filters, using common hne ratios can underestimate the hne flux, as the 
additional lines do not enter in the filter. 

5. We study the magnitude-color diagram used to select the candidates and we compute 
an analytic form of the standard deviation of the color of the objects. The standard 
deviation depends on a number of parameters, namely the exposure times in both 
bands, the colors of the object and of the background, the surface brightness of the sky 
and the magnitude of the object. The equation of (y[Q\ can help when there are not 
enough objects detected to calculate the standard deviation from the data available. 

6. We show how to compute the redshift range where a candidate can be included in 
the sample. This range is used to compute the comoving volume associated with each 
object. Using a global comoving volume for all the object tends to under estimate the 
volume for faint objects and objects selected near the selection curve. 

7. We provide equations to compute the optimal exposure time in the narrow band given 
the exposure time in the broad band. Given a color ratio for an object, we show the 
equation that relates the exposure times in both bands needed to effectively select the 
object. 

8. Additionally, if we impose the condition of minimum total exposure time, we obtain 
directly the exposure times in both bands. These equations are invaluable tools when 
planing narrow-band observations, as exposure time calculators usually only manage 
broad-band filters. 
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